5th generation 
machines emerge. 


exceeds $50. 


The pocket 

computer The pocket 
appears, computer 
priced at drops to 
$1000. $500 price. 


& 


Cost of random 
access storage 
falls below 


A new 

largest  °001f per bit. 

prime 

number 

is Cost per line of 
found. complex software 


exceeds $75. 


& The number of 
persons employed 
full time as 


"programmers" personal 

starts to drop. computers 
in U.S. 
exceeds 
250,000. 


Annual Forecast 
for the Coming 


Decade in Computing 


April 1976 Volume 4 Number 4 


Demise of PL/I 
becomes evident. 


Cost per line of 
complex software 


The‘ time-sharing 
industry dies. 
4 


Number of 


Big Brother arrives 
two years early 
through government 
surveillance via 
EFTS. 
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Number of makers 
of large main 
frames drops to 5. 


Computer control 
of automobiles 
begins. 


The pocket calculator 
industry hits 35 
million units per year. 
/ 


Cost of executed 
instructions falls 
below 100 million 
per dollar. 


More than half 
the computing 
power of the 


U.S. is in 
Price of the small dedicated 
pocket computer machines. 


drops to $250. 


The ABA forms 
Cost of executed . an active 
instructions falls 
below one billion 


per dollar. 


Computing courses are 
part of the curriculum 
in more than half of 
all U.S. high schools. 


Cost per line of 
complex software 
Cost of random exceeds $100. 
access storage 
falls below 


-000l¢ per bit. 


Cost of executed 
instructions falls 
below 10 billion 


Top speeds . 
of large per dollar 


machines 

reach on : 

Aa aecore Last traditional 
add time. textbook publisher 


gives up. 


computer group. 


This month's cover features our annual Forecast 
of events in computing for the coming decade. 
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The Forecast is prepared by the staff of POPULAR 
COMPUTING . it has been noted that long range forecasts 
tend to be pessimistic (that is, the predicted events 
come about sooner than anticipated) whereas short 

range forecasts are optimistic (that is, the time flies 
by and nothing happens). The crossover point seems 

to be around six years. 
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For example, 
the numbers on the right show the successive contents of a 
block of storage addressed at T, Trl, T+2, etc. 


Stack Action 


A main routine outputs a series of numbers that are 
integers in the range from 001 to 999; they are produced 
randomly. 


A subroutine is to place these numbers in a stack 
(such as a block of consecutive words of storage) in 
ascending order, with the number 1000 at the end of the 
stack as a terminator. The process will end when there 
are 50 numbers in the stack. 


There ig one restriction: no number shall be added 
to the stack if its value is within 10 of any previously 
stacked number. For example, if the number 383 has been 
placed in the stack, then any number in the ranges 
373--382 or 384--393 (and, of course, a repeat of 383) 
is to be rejected (i.e., by-passed). 


for the sequence produced in the left column, 
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start 1000 

383 383 1000 

572 383 _572 1000 

911 2283) 2572 9 911 1000 

100 D200 43880) = S72), 1912) 2008) 

582 reject 

250 100+ = 250, 383. . bya agnaeeiang 

908 reject 

110 reject 

398 100 _250 _383 _398 _572 _911 1000 
ees T+e TH3 THY T+5 TH+6 

What is the logic of the action called for by the 

subroutine? 
Outline a procedure to test that logic. LJ 
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Contest No. © 


In the article "Reading, Writing, and Computing" 
in the November-December issue of Creative Computing, 
Walter Koetke presents this problem: Is each positive 
integer 1 through N a divisor of some integer that 
contains only the digits one and zero? 


The answer is yes, as shown by David Ferguson: 


Any number N, not divisible by 2 or 5, divides 


N 
role. 4. If N = 225m, then N divides 


(10(™) =o), 2 yomax(asb) - MN). (x) is the Euler 


function: the number of numbers less than x relatively 
prime to x. 


A list of such numbers, for the integers from 1 to 
40, is given in the accompanying table. Any other digit 
could be substituted for the digit one. 


A possible algorithm for finding such numbers is 
given, using N = 29 as an example. We start with the 
four possible leading three digits (100, 101, 110, 111) 


and record the remainders on division by 29. Each of 
these remainders can be extended by appending a zero 
or a one, leading to a new set of remainders. The 


process continues as shown until a remainder of zero is 
found, after which the sequence leading to that zero 
(shown by the arrows) is readily traced. There is no 
assurance that this algorithm will give the smallest 
possible result. 


The suggested algorithm is probably not efficient, 
and would be awkward to implement for large values of N, 
since the list of remainders to be kept track of is of 
the order of N itself. 


For our 6th contest, the usual $25 prize goes to 
the person submitting the longest list (produced by 
computer, of course) like that of Table A. The list 
must be of consecutive integers, and must include the 
other factor, as shown. A flowchart (or equivalent) 
of the logic used to derive the results must be included. 
All entries must be received by POPULAR COMPUTING by 
June 30, 1976. 
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100 13 50 el 


101/14 3 si >|” 
= Ng 


180 
181 tf 
190 16 
191 17 37969 
29 | 1101101 
200 *— 26 
201 27 
220 ILf/ 
2el 18 
30 1 
31 2 
A possible algorithm 
4o ll for deriving the 
number composed of 
41 12 only ones and zeros 
that is divisible 
210 16 by 29. 
ali 17 
70 12 
71 13 
160 15 
161 16 
80 22 170 = 25 
81 23 171 26 
90 3 260 28 
91 4 261 ie) . 


PAINTED BY 


paint ESAMOn ey, 


Aigeld 


17002 Venture Bilvd., Encino, Ca. 91216 
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N times 
1 1 
2 5 
3 37 
4 25 
5 2 
6 185 
7 1443 
8 125 
9 12345679 
10 1 
11 1 
12 925 
13 (ai 
14 715 
15 74 
16 625 
ae 653 
18 61728395 
19 5209 
20 
21 481 
22 
23 43957 
a4 4625 
25 4 
26 385 
27 411522593 
28 3575 
2g 37969 
30 
31 322581 
32 3125 
33 336397 
34 3265 
35 2 
36 280558641975 
37 
38 2895 
39 259 
ko 25 


product 


1000 


11101 
1111111110 
100111 

100 

10101 

110 

1011011 
111000 

100 

10010 
11111110011 
100100 
1101101 
1110 
10000011 
100000 
11101101 
111010 
10010 : 
10100111111100 
111 

110010 
10101 

1000 


Table A. The first 40 integers, 


with numbers they divide, these 


numbers made up solely of i's and O's. 
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PROBLEM 1 2 5 


In essay No. 8 in the Art of Computing series a 
(PC28-3) the problem posed was this: / So/ 
@ Vi 
on 
A man cashes a check at a bank. The teller 
mistakenly interchanges the amount for dollars 
and the amount for cents, so that the man receives 
more than the amount of the check. The man 
finds that, after spending $14.19, he still has 
twice as much as the check cailed for. What 
was the amount on the check? 
In attempting to solve the problem algebraically, t 
we were led to the equation = 
oO 
A 


199D - 98C + 1419 = 0 (A) 


Since this single equation has two variables, ordinary 
methods of algebra will not do for solving it. If we 
solve the equation for C, we have 


¢ = 199 p + i419 
98 98 


G which can be written (separating out whole numbers): 
4 
C=) 2D +36 D+ 14 oe (B) 


The second condition on equation A (and all subsequent 
equations) is that each variable can take on only integral 


values. In equation B, then, the indicated terms are 
intrinsically integral, and hence the other two terms must 
also represent integers. This fact can be represented as: 
3 47 
== Ds Se 
98 98 
or 3D + 47 = 98E 
and, solving for D: i) 8 98 jn AT 
3 3 
Separating out the known integers as before: 
2 2 
J) S25, sp Sa NG) 
apa ig 


and this process can be repeated: 
2 -2 


€ =8 sates : 


i) 
Bi 
' 
wo 
it 
uy 
ny 
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until we reach a stage at which suitable integer values 
can be found. In this case, F = 0 will do nicely and, 
substituting back, we find that D = 17 and C = 49, 
The process just described is awkward to apply to 
the general case with literal coefficients, as in the 
third problem presented in issue No. 28: 


A man cashes a check at a bank. The teller 
mistakenly interchanges the amount for dollars 
and the amount for cents, so that the man receives 
more than the amount of the check. The man 


finds that, after spending X dollars, he still 
has twice as much as the check called for. What 
are the values of X that make this a real problem? 


However, we can go over the specific case and gain some 
insights into a general solution. The following material 
comes from Jules Mersel. 


Equation (A) may be expressed as follows: 
199D + 1419 = O mod 98 


3D + 1419 = O mod 98 


The first insight comes from noting that 1419 is a 
cumbersome number which may be reduced to its lowest 
possible value modulo 98: 


3D + 14(98) + 47 =O mod 98 
3D + 47 = O mod 98 


Similarly, since 47 is not a multiple of 3, we may force 
the relation: 
3D - 51 = O mod 98 
D = 17 mod 98 
and, since 0O=D<99, D = 17. 


Working Exes now, with the value 17 for D, 


we find that C = 49 and the original check (as before) 
was $17.49. 


In general, then, we have 


1000 + D - X = 2(100D + Cc) (1) 
199D + X = 98¢ (2) 
199D + X = O mod 98 (3) 

3D + X 2 O mod 98 (4) 


and our insight gained from the special case lets us 
write for X: 


X = 98m + r where O=r<97 


In turn, r=3n+p where 0= n= 32 
and 2 2 OO, ly Ge 2 
or X = 98m + 3n +p (5) 


where O=> m=102 
3D + 98m + 3n + p = O mod 98 
3D + 3n + p = O mod 98 (6) 
and there are three cases: 
Case I: p=l 


Case II: p = 2 
Case III: p = 0 


3D + 3n + 1» O moa 98 
3D + 3n + 99 = O mod 98 
D = -n - 33 mod 98 
D = 65 - n mod 98 
and since 0=D=99 and O=n = 32 
D=65-n 
Applying this value to equation (2): 
g8c = 199(65 - n) + x 
= 199(65 - n) + 98m + 3n +1 
= 98m - 196n + 12936 
= 98m - 196n + 98( 132) 


but O=C=99 therefore O=m -2n + 132299 
or 132=@n - m=33 


since n=32 and m=O, 
a solution for Case I exists only when 
an - m= 33 
and then D= 65 -n 
C = 132 + m = en 
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For example, choose n = 20, m = 4: 
x = 4(98) + 3(20) + 1 = 453 @ 
D = 65 = 20 = 45 
C = 132 +4 - 40 = 96 

and we have 96.45 - 4.53 = 91.92 = 2(45.96) 


Case II (p = 2): 
From equation (6): 
3D + 3n + 2 = O mod 98 
3D + 3n - 96 = O mod 98 
D+n - 32 = 0 mod 98 
D # 32 - n mod 98 
since O=[ D=99 and n=O, 
D= 32-n 


and, applying this value to equation (2): 


98c = 199D + x 
= 199(32 - n) + 98m + 3n + 2 
= 98m - 196n + 6370 
C =m - 2n + 65 


but O= C=99 therefore O=m - 2n + 65=99 
or -652m - 2n= 34 
since n= 32, there is a solution to Case II only when 
Qn - m= 34 
D= 32 -n 
C = 65 +m - 2n 


For example, choose n = 30 and m 20: 


x = 20(98) + 30(3) + 2 = 2052 


De s2 56 sa) 4 2 @ 
C = 65 + 20 - 2(30) = 25 [ 


and 25.02 - 20.52 = 4.50 = 2(2.25) 


Case III (p = 0): 
From equation (6): 3D + 3n # 0 mod 98 


D # -n mod 98 since 0O=D=99 


III(a) D=Oifns0O 
III(b) D=98 -n if n=0 
Applying D = 0 to equation (2): 


98c = 199(0) + x 
= 199(0) + 98m 
C=m 
So a result exists for Case III(a) when m=99. 
Applying D = 98 - n to equation (2): 
98c = 199(98 - n) + 98m + 3n 
= 98m - 196n + 98(199) 
C =m - 2n + 199 


since 0=C=99, therefore O=m - 2n + 199599 


but, since 1=n=32 and O= m=102, there is no 
solution to Case III(b). 


So for Case III (p = 0), a solution exists only for 
n = O and m=99 


and then DB ©) 


As an example, choose m = 20: 
x = 98(20) + 0 + 0 = 1960 
D=0 
C = 20 
and 20.00 - 19.60 = 0.40 = 2(0.20) 
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A If you need to 
tt SOMEWHERE IN HERE... keep up with 
af the computer field e 
au there is an article but lack the 
about computers time to research 
that may be and read, 
helpful to you. 
DATA 
These are the 203 PROCESSING 
business and DIGEST 
computer journals 
we read is your answer. 
every month... Monthiy, averaging 
and digest the best 50 items, including 
of the articles book reviews; index, 
for calendar, complete 
YOU references to 


original articles. 


Write for information to Dept D, 


Data Processing Digest,Inc. 3% 


6820 LA TUERA BOULEVARD. LOS ANGELES, CALIFORNIA 90045 / PHONE (213) 776-4334 


Jared Weinberger, Bologna, Italy, points out that the 
sequence given as an illustration for SLOWGROW (PC30-10) 
is in error. The sequence starts with any integer and 
forms each subsequent term by adding to the preceding term 
the sum of its digits. Thus, the example given should 
have begun 


7 1 19 29 40 44 52 59 73... 


Log 37 1.5682017240669949968084506895 3912944798297 2690166313 
In 37 3.610917912644224444 36809567 1031447 16 3900077 587 167636 
V37 6.082762530298219688999684.24520206706208497 0094786411 
V3T 3.332221851645953260095450505185119004409616671950062 
V37 1.43489516567535945282940584.2182258706024214075867060 
‘V3T 1.036769033849400879875881175233017497034970001106672 


es! 11719142372802611 . 30877293979119019452167536 369446182 
23834809114 359096826486861 321089074 


2480534602660760780 .433710790045522100100824031419549 
833903269 349676966827477951 312798 


tan ~ 37 1.5437758776076 318 3044 3146 3582812318099037 33691908670 


N-SERIES yf 


e (MO[N/O/P/o}L| | |? 


...DWalderdas ! 


IT'S TIME TO TRY Uy, Your Word 


e THE NEWEST, MOST CHALLENGING WORD GAME 


For two or more players. 


ov 
Bj will both challenge and entertain you and your friends. 
3) | Increase your word skill--organize a tournament! 


:0) 


> 
"Slo 


The complete game, along with instructions 
and several variations, is only $2.00! 


ORDER FROM Popular Computing 


BOX 272 CALABASAS, CA. 91302 


Ey 


iw 
=p] 


Up Your Word — Conrtight 1976 by Poplar 


l For each game, send $2 plus 30¢ for postage and handling 
1! U.S., Canada, and Mexico) 
7 or plus 50¢ for postage and handling 
all other countries) 
For two games, send $4 plus 50¢ (U.S., Canada, and Mexico) 
or plus 75¢ (all other countries) 


California residents add 12¢ per game. 
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Life or Death 


The following game is to be played. A number is 
entered in one cell of an infinite grid. At each stage 
of the game, the number in a cell dictates a move in one 
of 8 directions: 


1 2 3 

4 x 5 

6.47 io 
according to X modulo 8. For example, if the number in 
cell X is 5, or 13, or 21, or 69, the move is made to the 
east. To make the move, the value of X is reduced by 
one, and the cell it points to will contain either X/2 
(if X was even) or 3X+1 (if X was odd). No further 


movement arises from a cell containing the value one. 


All moves for a stage are made simultaneously. When 
two or more cells are all creating new numbers for the same 
cell, the lowest number prevails. 


if the game begins with the number one, it is 
immediately over. If the game begins with 2, it is over 
in one move. Things get interesting when the game begins 
with 3. The accompanying figure shows the first ten 
stages for the game beginning with 3. It can be seen 
that large blocks of ones develop, since it is well known 


that the 3X+1 process always converges to one (see the 


discussions of the 3X+1 process in issues 1, 4, 13, and 25). 


We have then these problems: 


(1) Write a program to carry out the Life or Death 
procedure from stage to stage. 


(2) Will any pattern stabilize; that is, will any 
pattern converge to all ones? 


(3) What 1s the developing pattern for starting 
values other than 3? 


St~LE9d 


Life or Death 
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90:89 -88 °87 -86-85°84-83-82-81-80 


Problem Solution 
The IBM Card Problem 


Problem 55 called for: 


A. The logic of a subroutine that furnishes the 
number of paths from hole (1,9) to hole 
(column X, row Y) of a standard IBM cara. 

B. A Fortran program for that logic. 

C. The number of paths from (1,9) to (80,12). 


Contributing editor Thomas R. Parkin writes: 


NUMBERSS5 


Let us discuss this problem in terms of the vertical 


and horizontal segments of the paths from one point to 
another. All paths must proceed from left to right and 
from bottom to top, only. Thus we may assign O to each 
vertical segment from one row to another and 1 to each 
horizontal segment from one column to another. Thus 

to go from (1,9) to (3,8), we can go across 2 columns 


and up 1 row, or across 1 column and up 1 row and across 


1 column, or go up 1 row and across 2 columns--for a total 
of 3 possible paths. These paths correspond to 110, 101, 


and O11, respectively, according to our labelling. Thus 
we see that if we go from column A to column B and from 
row C to row D, we will have (B - A) 1's and (D - C) O's 
in our final path descriptions. (Of course, we have 
converted the designations of rows on an IBM card from 
Cpe Ole Oen, s,9) Into de, 15 10, (O02 .%, 2) 


Let us say we traverse N columns (N = B - A) and 
K rows (K = D - C). Then we have exactly as many paths 
as we can form binary numbers each consisting of exactly 
N 1's and K O's and, of course, this is 


N (N+K) ! 
P Sa = 
(N,K) N+K N! K! 


Thus, if we wish to proceed from (1,9) to (3,8), we are 
going up 1 row and across 2 columns, so N = 2, K = 1 and 


P(n,K) S P(2,1) = 3!/1!2! = 3 


as we found earlier. Similarly, to go from (1,9) to 
(80,12), we have N = 79, K = 11, and 


P(79,11) = 90!/79!11! = 


= 41604694413840 
1°2°3°4°5°6°7°8°9-10°11 


The pertinent Fortran program follows. 


AAADQAARQAARARARQARAARAAAAAAAAAA 


AARAQANARQRIaA0a 


aaa 


aQAAandgana 


SUBROUTINE PATHS (IC, IR, IP, AP, ILM) 
This subroutine calculates the number of paths 
from hole to hole in an 80 x 12 hole card, 
from bottom row, left column (or IC = 1, 
IR = 9) to the hole at (IC, IR), expressed 
in standard IBM card numbering conventions. 
The paths are constrained to go from left to 
right and bottom to top only. Columns are 
numbered 1 to 80 left to right, rows are 
numbered 12, 11, 0, 1, 2,...,9 from top to 
bottom. 


Calling sequence: CALL PATHS (IC, IR, IP, AP, ILM) 


where IC is the column number and 
IR is the row number of the terminating 
point. ; 
IP will contain the number of paths between hole 
(1,9) and hole (IC,IR) upon return, provided 
the integer limit of the particular Fortran 
being used can accommodate the number; if not, 
IP = -l, and the floating point value of the 
number of paths will be in AP. 
ILM is preset to the (Fortran integer limit -1). 


The following statements are inserted to skip 
the calculations if parameters outside the 
proper limits are given. 


IF (IR.EQ.10) G@ TS 50 

IF (IR.LT.00) G# TP 50 

IF (IR.GT.12) Gd TS 50 

IF Weereecaa) Gg TS 50 

IF (IC.GT.80) GZ TS 50 
These could be combined into one compound 
logical statement if the particular Fortran 
will allow this form. 
We will now convert the input parameters to 
internally useful ranges. We shall let K run 
from 0 to 11 corresponding to inputs of rows of 
9,8,...,0,11,12, and N run from 0 to 79 
corresponding to columns of 1 to 80. 

K = 9 - IR 

IF (IR.GE.11) K = IR-1 

N=IC-1 
Now check for special cases of N or K or both 
being zero. 


IF eee GO TS 30 
IF (K.EQ.0) GG TS 40 


If we fall through to here, neither N nor K is 


zero, and we can proceed. In order to calculate 


the number of paths, we must evaluate the 
formula P{N,K) = (N+K) FACTORIAL DIVIDED BY 
N FACTORIAL TIMES K FACTORIAL. 
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AaAagaaaaga 


ARAAARQRAKDRAARAAaAaAQa 


AQAAAa 


10 


20 


30 
4O 


50 


In order to shorten the calculation as much as 


possible, we will first find out which is larger, 


K or N, and cancel off that many factors from 
both the numerator and denominator, 

To do this, we will preset the limit of the 
D@ LOOP index to one of the variables and 
reverse this if appropriate. 


IX = N 
IF (N.@T.K) IX = K 


The D@ LOOP limit is now the smaller of N or K, 


Now perform the evaluation of the factorial 


formula. Note that we have the same number of 


factors in both the numerator and denominator 
after cancelling the low order terms of the 

numerator against either N or K (whichever is 
larger) terms of the denominator. Furthermore, 
we shall work in floating point (we assume the 


floating point limit is greater than 4.2*10**12), 


We shall use a floating value of K or N 
(whichever is smaller) and our floating point 
answer space. 


IY = K 
IF (N.GT.K) IY = N 
AY = FLOAT (IY) 
AI = 0.0 
AP = 1.0 
DS 10 I = 1, IX 
AI = AI + 1.0 
AP = AP¥(AY + AI)/AI 
CONTINUE 
Our answer is now in AP. We must check the 


integer limit and tidy up the special cases 
of N or K or both equal to zero. 


B = FLOAT (ILM) 

IF (AP.GT.B) G@ TP 20 
IP = IFIX (AP) 
RETURN 

CONTINUE 

IP = =1 

RETURN 

CONTINUE 

IF (K.EQ.0) G@ TH 50 
CONTINUE 

IP 2l 

RETURN 

CONTINUE 

LPR s=" 0 

RETURN 

END 


